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SUMMARY 


The  present  paper  derives  Von  Karmmi's  Foiirler  Integral  method  in 
supersonic  wing  theory  directly  from  the  basic  concepts  of  the  harmonic 
source  and  doublet.  The  method  is  first  applied  to  investigate  the 
general  solution  of  the  wave  drag  of  a  tapered  swept  wing  with  a  sym¬ 
metrical  diamond  airfoil  profile.  The  general  solution  includes  all 
kinds  of  wing  plan  forms  which  may  be  swept  backward  or  forward,  and 
tapered  or  reversely  tapered  to  any  ratio.  A  number  of  the  limiting 
cases  are  also  Investigated.  For  practical  aerodynamic  design,  two 
families  of  wing  plan  forms  with  the  fixed  taper  ratios  0.2  and  0.5,  any 
swept  angle,  aspect  ratio,  and  Mach  number  are  shown  in  graphs.  Some 
particular  applications  are  illustrated. 

The  reversed- flow  theorem  on  wave  drag  as  shown  by  Von  Karman  and 
Hayes  checks  well  with  the  consequence  of  the  general  solution.  This 
method  shows  a  certain  elegance  as  no  conical-flow  assumption  is  needed 
and  the  mathematics  is  powerful  enough  to  obtain  a  general  solution 
covering  all  possible  geometrical  arrangements  without  detailed 
considerations . 

While  in  recent  years  the  direct  problem  of  finding  the  lift  dis¬ 
tribution  with  a  given  angle  of  attack  of  the  wing  has  been  well  solved 
method  of  conical  flow  and  by  other  methods,  the  present  treat¬ 
ment,  on  the  other  hand,  investigates  the  inverse  problem,  that  is,  to 
find  the  downwash  distribution  in  the  plane  of  a  wing  with  a  preassigned 
pressure  distribution.  This  is  particularly  favorable  with  the  present 
method.  The  general  solution  of  the  downwash  of  the  tapered  swept  wings 
is  derived  for  the  case  where  a  constant  lift  distribution  on  the  wing 
is  preassigned.  Of  course,  the  method  may  be  applied  to  any  lift  or 
pressure  distribution  along  the  wing  chord  and  span.  The  corresponding 

of  attack  of  the  wing  and  the  downwash  can  be  determined  everywhere 
in  the  plane  of  the  wing.  To  demonstrate  the  downwash  distribution  as 
given  by  the  general  solution,  the  downwash  of  a  number  of  wings, 
including  a  sweptback  tapered  wing  with  supersonic  trailing  edge  and  a 
triangular  wing,  is  shown  in  graphs. 
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The  general  solution  of  the  direct  problem  in  the  present  formula¬ 
tion  is  a  singular  integral  equation  of  the  Weiner-Hopf  type  and  is 
difficult  to  obtain. 


INTRODUCTION 


In  1935  Von  Karman  first  shoved  that  the  concept  of  the  Fourier 
integral  can  be  adopted  to  explain  the  similarity  of  Prandtl  ving  ^theory 
and  the  theory  of  planing  surfaces  (reference  l).  In  19^6^  Von  Karman 
introduced  the  Fourier  integral  method  to  the  supersonic  wing  theory 
(an  unpublished  report  of  Northrop  Aircraft^  Inc.)*  The  present  author 
under  Von  Karman*  s  guidance  investigated  this  problem  to  a  certain  extent j 
some  of  the  results  have  been  published  (reference  2)  but  some^  presented 
in  19^  at  the  Sixth  International  Congress  for  Applied  Mechanics,  are 
not  available  in  published  form.  The  present  report  may  be  considered 
as  an  extension  of  the  earlier  work.  As  some  of  the  literature  is 
unavailable  at  this  moment,  some  of  Von  Karmen’s  work  is  repeated. 

As  far  as  the  linearized  theory  on  supersonic  wings  is  concerned, 
the  conical-flow  method  was  originated  by  Busemann  (reference  3) 
later  reproduced  and  extended  by  Lagerstrom  (reference  k)  ^  Hayes 
(reference  5)^  Stewart  (reference  6),  Laporte  and  Bartels  (reference  7); 
and  Snow  (reference  8).  Each  of  the  later  investigators  approached  the 
same  method  with  slightly  different  techniques  but  confirmed  the  essential 
results  of  Busemann.  Before  the  conical-flow  method  was  known  in  this 
country,  Jones  visualized  the  advantages  of  the  conical  flow  and  showed 
some  of  the  basic  physical  concepts  (reference  9)  1^^  19^5*  Later,  in 
191^-6,.  with  the  concept  of  conical  flow,  Jones  also  showed  the  invariance 
of  the  Lorentz  transfomation  and  introduced  the  oblique  coordinate  in 
the  swept-wing  problem  (reference  lO) . 

The  challenge  of  supersonic  flight  has  aroused  the  interest  of  many 
other  Investigators.  Puckett  showed  that  SQurce  integration  method 
could  be  applied  to  study  the  wave  drag  problem  (reference  11).  Later 
Eward  extended  this  method  to  solve  the  lift  problem  of  a  finite  wing 
of  any  plan  form  (references  12  and  13)*  Heaslet,  Lomax,  and  Jones 
extended  Volterra'  s  and  Hadamard's  method  to  the  supersonic  wing  problem 
(references  ik  and  I5) .  Gunn  applied  the  operational  calculus  to  the 
same  problem  (reference  16) .  There  are  many  other  works,  all  of  which 
cannot  be  mentioned  in  this  paper. 

On  the  other  hand,  the  basic  concept  of  the  Fourier  integral  method 
is  quite  different  from  the  above  methods.  Instead  of  the  concept  of 
conical  flow  or  simple  sources  and  doublets,  the  present  method  considers 
along  the  direction  of  flight  infinitely  long  hamonic  source  (or  doublet) 
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lines,  the  behavior  of  which  is  quite  equivalent  to  the  harmonic  acoustic 
source  (or  doublet)  in  the  sequence  of  time.  Each  of  such  source  lines 
will  send  out  a  divergent  cylindrical  wave  in  the  radial  direction.  The 
potential  of  such  a  wave  can  be  expressed  by  the  product  of  the  source- 
strength  function  and  the  Hankel  function  of  the  second  kind.  For  each 
eigenvalue  or  frequency  of  the  source  oscillation,  there  is  a  corre¬ 
sponding  eigenvalue  in  the  argument  of  the  Hankel  fiinction.  By  means 
of  the  principle  of  superposition,  an  arbitrary  distribution  can  be 
synthesized  with  such  simple  harmonic  sources  of  different  frequencies, 
if  the  frequency  spectrum  is  a  continuous  one.  The  most  powerful  tech¬ 
nique  to  serve  such  a  purpose  is  the  Fourier  integral  method. 

Now  in  the  finite-wing  problem,  the  boundary  is  considered  as  a 
source  sheet  rather  than  as  source  lines.  Such  a  sheet  can  be  built  up 
by  integrating  elements  of  a  source  strip  which  is  equivalent  to  a 
source  line  in  behavior. 

The  above  simple  physical  interpretation  may  help  to  give  some 
insight  into  the  mathematical  theory  to  be  discussed  later. 

The  present  paper  first  tries  to  give  an  introduction  to  Von  Karman's 
method  in  terms  of  the  physical  rather  than  mathematical  approach.  Then, 
the  wave  drag  of  a  tapered  sweptback  wing  is  treated.  This  work  is  not 
a  duplication  of  that  of  some  other  concurrent  investigators,  but  rather 
a  supplementary  contribution.  The  next  step  is  the  lift  problem.  For  a 
given  angle  of  attack,  on  an  arbitrary  wing  plan  form,  it  is  rather 
difficult  to  solve  the  integral  equation  analytically.  But  for  a  given 
lift  distribution  the  downwash  angle  can  be  evaluated  anywhere  in  the 
plane  of  the  wing.^ 

The  method  has  a  certain  elegance.  The  complete  physical  effect 
of  the  arbitrary  wing  can  be  expressed  in  terms  of  one  integral.  The 
finite  number  of  discontinuous  points  of  the  Fourier-Bessel  integral 
gives  exactly  the  right  picttire. 

This  investigation  was  conducted  at  the  Department  of  Aeronautics 
of  the  Johns  Hopkins  University  under  the  sponsorship  and  with  the 
financial  assistance  of  the  National  Advisory  Committee  for  Aeronautics. 
Due  appreciation  should  be  given  to  Misses  V.  O'Brien  and  P.  darken  for 
their  assistance  in  carrying  out  the  research. 


^In  some  recent  unpublished  Navy  reports,  J.  D,  Miles  has  obtained 
a  solution  of  the  direct  problem,  particularly  for  the  rectangular  wing 
tip  with  constant  angle  of  attack. 
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GENERAL  THEORY  OF  THE  FOURIER  INTEGRAL  METHOD 
Elementary  Solutions 


In  the  three-dimensional,  steady  supersonic  flow  of  a  compressible 
nonviscous  fluid,  the  differential  equation  of  motion  in  the  linearized 
sense  Is 

0yy  +  0ZZ  = 


where  the  velocity  potential  0  is  defined  by 

grad  0  =  q  -  Ui  =  ul  +  v3  +  wk,  the  disturbance  velocity  vector.  It^ 
may  be  considered  as  a  two-dimensional  wave  equation  if  x  is  conceived 
as  time  in  the  sequence  of  which  the  futiire  can  contribute  nothing  to 

the  present.  This  can  be  shown  clearly  by  Introducing  t  = 
that  is^ 

0yy  +  0ZZ  =  ht 

With  the  dimension  in  the  flow  direction  being  equivalent  to  time  t. 
After  this  transformation,  the  velocity  of  propagation  in  equation  (2)  is 
unity.  In  this  sense,  equation  (2)  becomes  the  potential  problem  of  the 
acoustic  source  and  doublet  in  the  y, z-plane  in  the  sequence  of  time  as 
pointed  out  by  Von  Karman  (reference  2).  After  this  transformation,  the 
plan  form  has  to  be  ineadjusted  so  that  the  Mach  waves  are  inclined 
backward  at  45^  from  the  flow  direction.  As  an  example  of  transforma-  . 
tion,  figure  1  gives  a  sweptback  wing  in  the  physical  and  transformed 
planes . 

For  a  simple  time -dependent  harmonic  source  with  strength  cos  Xt 
(where  x/2jt  is  the  frequency  of  oscillation)  located  at  (ti^O)  in  the 
y, z-plane,  the  elementary  solution  of  equation  (l)  as  given  by  Lamb 
(reference  17^  P-  297)  is 

0o(t,y,z;  11,0)  =  Rej^-  I 


where  r  =  -  T))^  +  z^ 

second  kind  of  zero  order. 


and  HQ^^^(Xr)  is  the  Hankel  function  of  the 
It  is  easy  to  show  that  the  solution  satisfies 


MCA  TN  2317 


5 


both  the  differential  equation  (2)  and  the 
iX.t 

cos  Xt  =  Re  e  .  At  very  large  values  of 

0o(t,r)  s;  ^  cosrx(t 

^8jtXr  L 

the  amplitude  of  which  is  inversely  proportional  to  the  square  root 
of  r.  Physically^  this  represents  the  potential  at  the  point  (y^z)  at 
time  t  due  to  a  harmonic  source  of  strength  cos  \t  at  distance  r 
from  it .  Such  a  disturbance  potential  is  called  a  divergent  wave .  This 
disturbance  propagates  from  the  point  (t],0)  with  a  circular  or  cylin¬ 
drical  symmetry. 

Now  if  the  strength  of  the  source  is  the  real  part  of  F(X)ei^*^ 
instead  of  just  e^^^,  the  potential  must  be  also  written  as 

-^,0)  =  i  F(X)e^^^Ho^^^(Xr)]  (4) 

where 

F(X)  =  Fq(X)  +  iFi(X)  (5) 


given  boundary  condition 
Xr  equation  (3)  becomes 


with  Fq(-X)  =  Fq(X)  being  even  and  Fi(-X)  =  -F]^(X)  being  odd.  It 

is  understood  that  F(\)  is  a  function  representing  the  strength^ 
amplitude^  and  location  of  the  source  line. 


If  two  such  harmonic  line  sources  of  equal  but  opposite  strength  - 
the  negative  one  located  at  (t],-^/2)  and  the  positive  one  at  (t],^/2)  - 
approach  each  other^  a  hamonic  doublet  line  can  be  obtained  if 
lim  F(X)^  =  G(X)  is  considered  as  a  finite  quantity: 

0l(t,y,z}  Ti,0)  =  0Q^t,y,zj  r\, 


=  llzn 


=  Re 


(6) 


6 


MCA  TN  2317 


where  the  doublet  Is  defined  as  positive,  and 

G(X)  =  Gq{X)  +  IGi(x)  (7) 

with  Gq(-X)  =  Gq(X)  being  even  and  G^^C-X)  =  -G2(X)  being  odd.  It  is 

very  interesting  to  see  that  this  potential  of  the  doublet  no  longer  has 
cylindrical  symmetry^  but  is  anti  symmetrical  vith  respect  to  the 
t^y-plane.  The  harmonic  doublet  line  will  be  used  for  the  wing  lift 
problem. 


Boundary  Conditions  in  the  Supersonic  Wing  Problem 

Let  a  flat  body  or  wing  occupy  a  region  in  the  x,y-plane  and 
and  a_  be  the  slopes  of  the  upper  and  lower  wing  surfaces  at  (x^y), 
respectively.  Both  are  small  in  comparison  with  unity  and  are  zero 
outside  the  finite  region  occupied  by  the  wing.  Within  the  approxima¬ 
tion  of  the  linearized  theory^  the  wing  may  be  considered  equivalent  to 
the  superposition  of  a  symmetrical  body  and  a  mean-cambered  surface. 

The  slope  of  the  symmetrical  body  is 

“o  =  |(a+  -  a_)  (8) 

and  the  slope  of  the  cambered  surface  is 

ttp  =  |(a+  +  a_)  (9) 


where  and  are  considered  at  the  x,y-plane.  In  other  words, 

the  effect  of  exact  location  of  and  in  the  z-direction  is 

entirely  neglected  in  the  future  treatment  which  has  been  shown  by 
Von  Kann^  to  be  consistent  with  the  linearized  theory  for  the  flat  body. 
From  the  above,  the  required  boundary  condition  of  the  potential  is  the 
normal  velocity  to  the  wing  surface  which  is  zeroj  that  is, 

qn  =  0  (10) 

In  this  case  of  a  flat  body,  where  the  vertical  velocity  conrponent  pre¬ 
dominates,  equation  (lO)  can  be  approximately  expressed  by 


(11) 
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where  w+  and  w_  represent  vertical  velocity  components  on  the  top 

and  bottom  of  wing  surfaces,  respectively.  Besides,  at  a  distance  far 
away  from  the  wing  in  the  upper  stream,  u=v=w=0  if  there  is  no 
other  distiorhance  generated  in  the  upper  stream. 


Supersonic  Flow  about  a  Flat  Symmetrical  Body 

If  a  source  line  at  (ti,0)  is  of  any  arbitrary  strength  with  finite 
time  t  the  potential  can  be  built  up  with  equation  (!<■)  by  means  of  a 
Fourier  integral  as 

Hoo 

=  ReJ  dx|^-  F(X,Ti)e^^''^Ho^^^(Xr)j  (13) 

where  U  is  the  free-stream  velocity.  The  insertion  of  U  in  the 
source-strength  fvmction  is  purely  for  future  convenience. 

The  above  Integral  and  its  first  derivatives  exist  if  the  Lebesgue 
integral 


X  |F(X,T))Ho^^^(Xr)l  dX.  <  00 


The  above  condition  is  automatically  satisfied,  if  the  soinrce  with  finite 
strength  is  in  action  for  a  finite  time  interval.  Actually,  equation  (13) 
can  be  considered  as  the  potential  of  a  body  of  revolution.  Sometimes, 
in  conforming  with  the  convention  of  the  complex  Fourier  integral,  the 
above  equation  may  be  conveniently  written  as 


0o'(t,y,z)  =  -  Y  dXF(X,Ti)e^^%0^2)(^^) 


(1^) 
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with  the  understanding  that  the  path  of  integration  of  X  is  slightly 
above  or  below  the  origin  in  the  con^jlex  \-plane  so  that  the  singularity 
at  X  =  0  can  be  detoured.  It  can  be  verified  easily  because 


F(-X,ti)  =  Fq(X,ti)  -  iF3^(X,Ti)  = 

Ho^^^(-Xr)  =  Ho^^^(xre^’')  =  Ho^^^(xre‘^’')  =  Ho^^^(Xr) 


where  the  bar  represents  the  complex  conjiigate  of  the  original  function. 
Now,  the  above  source  line  may  be  considered  as  an  elementary  strip  of 
a  source  sheet  in  the  t,y-plane  of  width  dr\  with  the  above  strength  as 
the  strength  density  per  unit  area.  It  is  clear  then  that  the  potential 
due  to  a  source  sheet  built  up  by  such  strips  is 


=  ;zSq*  dll 


•oo  noo 


=  _  dT)  I  dXF(X,Ti)e^^'*^Ho^^^(Xr) 


(15a) 


'  —  00  U-00 


For  the  convenience  of  the  later  development,  f(X,T))  =  -  ^  F(X, t))  is 
introduced  into  equation  (l5a),  and  then 


»00  f 100 


0o(t,y,z)=^  dn  I  dXf(X,n)e^^'^Ho(2)(Xr) 


(15b) 


One  thing  should  be  pointed  out.  Owing  to  the  structure  of  the  integral, 
the  source  sheet  may  be  composed  of  a  number  of  discontinuous  portions 
in  the  y-direction,  and  the  source- strength  function  f(X,T])  may  be 
different  in  different  discontinuous  portions,  but  it  must  be  the  same 
within  one  continuous  portion  of  the  source  sheet.  Of  course,  in  that 
case,  the  integration  with  respect  to  r\  has  to  be  broken  up  accordingly 
Physically,  this  Integral  can  give  the  potential  not  only  due  to  one  wing 
but  also  to  a  number  of  wings  or  flat  bodies  in  the  t,y-plane.  There¬ 
fore,  it  can  be  used  to  study  the  interaction  between  wings. 


2. 
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'  The  preceding  equation  will  be  the  disturbance  potential  of  the  ving^ 
if  the  relation  of  the  source -strength  function  Re  can  be 

identified  with  clq  or  the  thickness  distribution  of  the  symmetrical 

flat  body .  •  This  can  be  obtained  from  the  boundaiy  condition.  Differ¬ 
entiating  equation  (l5b)  with  respect  to  there  results 


hz 


Ui 

2 


poo 


-00 


dT] 


X 


U-oo 


dXf  (X.r)  ^ 


(16) 


Now  as  z — ^0  the  integrand  goes  to  zero  except  in  the  neighborhood  of 

T)  =  y .  Let  T)  =  y  +  z  tan  0}  dq  =  z  sec^0  where  9  =  tan"*^  S-.-L  V  ^ 

z 


y  T] 


Besides,  «  2i/jtXr  when  X.r  or  r— yO.  Thus  equation  (l6) 

becomes 


(w)z=+0  = 


.Sz  j 


z=+0 


=  lim 
z — >+0 


Ui 


dXe^^^  I  X,f(X.^y  +  z  tan  9)  —  —  sec^0  d0 


\/2 


-rt/2 


TtX  p2 


U 


=  llm  ^  I  dXe^^'^  I  Xf{X,y  +  z  tan  e)de 


V/2 


z— >+0 


U 
jr 

u- 


«/2 

Jt/2 

dA,e^^‘^f(X,y)  j  d0 

-it/2 


=  U  I  dXe^^'^f(X,y) 


(17) 


-  00 
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From  equation  (lO), 


(^)z=+0  =  U(ao)z=+o  =  U  dXei^^^f(X,y) 


ao(t,y,+0)  = 


J|fo(X,y)  cos  \t  +  fj(X,y)  sin  X'^dX 
0 


From  the  theory  of  complex  Fourier  transfonn^ 


^  a(t,y)e"^^^  dt 


This  means  that,  if  the  source  distribution  f(X,y)  is  chosen  to  be 
the  Fourier  transform  of  the  distribution  of  the  angle  of  attack,  the 
potential  0q  will  represent  the  flow  about  the  wing  automatically. 

Since  the  angle  of  attack  is  zero  ahead  of  the  wing,  it  can  be  shown 
easily  that  equation  (l5b)  satisfies  the  second  boundary  condition  of 
equation  (12) . 


Supersonic  Flow  about  a  Lifting  Surface 

By  means  of  the  same  argument  as  in  the  last  section,  the  potential 
can  be  formulated  for  a  doublet  sheet  of  strength  Re  e^^'^G(X,T))  from 
equation  (6) : 

poo  poo 

0l(t,y,z)  =Re  dn  dx  X  G(X,Ti)ei^%l(2)(^j.)£  (20) 

J-  00  Jo 


where  a  constant,  free-stream  velocity  U  is  Introduced.  In  order  to 
simplify  the  later  development,  further  introduce 
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g(X,Tl)  =  ^  ^  =  go(X.,Tl)  -  igl(A,,Tl) 

where^  with  equation  (7)> 

go(^^’l)  =  — (even)  1 


gi(^^n)  =  — - 

t 

Then  equation  (20)  can  he  written 


( even) 

8m* 

XGq(X^  q) 

8m* 

(odd) 

Poo  p  CX) 

0l(t,y,z)  =  Re  -2m'U  dXg(X^r))e^^'^Hj_^^\xr)^ 

^  U-  00  (Jo 


The  above  equation  can  be  written  entirely  in  the  coii5>lex  form 

p  00  poo 

0l(t,y,z)  =  -m'U  I  dT)  I  dXg(>.,Ti)e^^%i^^)(Xr)^  (24) 


-  00  U"  00 


if  the  path  of  X  is  chosen  below  X  =  0  in  the  complex  X~plane^ 
because 


g(-x,il)  =  go(A.>Ti)  +  igi(X,^Ti)  =  g(X.,Ti) 


The  path  of  X  with  Re  X  >  0  cannot  be  chosen  in  order  to 
detour  X  =  0.  In  that  case^ 


Hi^^^(-Xr)  =  Hi(^)(xre^^)  =  -Hj^^^Xr)  -  2J2(Xr)  ^  H^(^)(Xr) 


which  is  impossible  to  reduce  back  to  the  original  equation  (23) 
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In  order  to  find  out  whether  equation  (24)  for  the  doublet  sheet 
satisfies  all  the  boundary  conditions  of  the  wing  or  not,  first  differ¬ 
entiate  it  with  respect  to  t. 


ht  ' 


dXg  ( X,,  T] )  e  ^  ^  ^  )  f 


(26) 


The  above  equation  is  identical  in  form  with  equation  (l6)  if 


f(X,‘n)  is  replaced  by  ^  g(^^T]).  Therefore,  with  the  same  technique, 
it  can  be  shown  that,  as  z — >+0,  T) — >y. 


/ z=+o  ^ 


=  ^  dXel«g(X,y) 


(27) 


Row  in  the  linearized  theory  the  pressure  coefficient  on  the  wing 
surface  is 


-2u  -2  r^i 


Cp(t,y,+0)  =  V  ~ 


’6jz5. 


U  Urn'  \at  /  z=+o 


(28) 


Thus,  with  equation  (27), 


Cp(t,y,+0)  =  dX,e^^'^g(X,,y) 


(29) 


With  the  Fourier  transform. 


g(X,y)  =  ^  I  dte^^^C^(t,y,+0) 


(30) 


Therefore,  if  the  pressure  distribution  along  the  wing  chord  is 
given,  the  doublet  distribution  fimctlon  g(X,y)  of  frequency  X/2n 
can  be  obtained  with  equation  (30).  The  function  g(X,y)  is  inde¬ 
pendent  of  the  pressure  distribution  of  the  neighboring  chord.  F\irther- 
more,  if  X->0  in  equation  (30)> 
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g(O^y)  =  ^  I  dtCp(t,y,+0) 


(31) 


On  the  right  side,  I  citCp(t,y,+0)  is  equal  to  one-half  of  the  lift- 


coefficient  distribution 


to 


"  2  V^y  J’ 


at  y.  Thus,  g(0,y)  is  equivalent 


at  y  which  has  a  physical  meaning 


Unfortunately,  equation  (23)  of  the  douhlet-sheet  potential  satisfies 
one  of  the  boundary  conditions  on  the  wing,  hut  does  not  satisfy  the  other 
required  boxmdary  condition  that  the  potential  and  its  derivatives  must 
be  zero  far  ahead  of  the  body.  In  other  words,  the  correct  potential 
must  be  zero  as  t— »-oo. 


The  potential  of  the  doublet  sheet  at  large  ±t  is  investigated  as 
follows.  Introduce  V  =  Xt  into  equation  (23): 


0(t,y,z)  =  Re 


>00  f  *00 


(32) 


For  very  large  values  of  t  there  can  be  written 

2it 


3tvr 


and 


Si  g(0,il) 


11m  =  lim  Re|  |  ^  dTig(0,ii) 


t--t±oo 


t  H^+oo 


=  lim 
t - ^±00 


-4m 'U1 

ij-oo 


.iv 


dv 


Ita’U 


dTig(o,n) 


£i2-Ydv 


7t 


0l(±“^y^z)  =  X*  =  ±2m'U  <iTig(0,Tj) 


(33) 
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where  the  negative  sign  corresponds  to  t  = • -00  and  the  positive  sign 
corresponds  to  t  =  00.  As  x*(y*z)  is  independent  of  t,  it  satisfies 

p.2  *  ;.2y* 

the  Laplace  equation  ^  ^  ^  =  0  and  can  he  added  to  0^  the 

doublet-sheet  potential^  without  an  effect  on. the  pressure  distribution. 
It  should  be  noted  that  the  value  of  the  potential  function  as  t — >00 
is  2X*.  Hence  may  be  interpreted  physically  as  the  potential 

function  of  the  downwash  in  the  Trefftz  plane. 

♦ 

Thus,  from  equations  (24)  and  (33)  the  complete  solution  of  the 
present  problem  on  supersonic  wing  lift  is 


0(t,y,z)  =  01  +  X*  =  m'U  dr\  2g(0,T])  -  I  r  dXg(X,,Ti)e^^%i^^^(X.r) 

.  r^ 


In  the  case  where  Cp  is  given^  the  above  equation  can  be  written, 
with  equations  (30).  and  (31)^ 


0(t,y,z)  =1^  -%dT)  dTCp(T,T],+0)  2  -  r  (Xr) 


J-co  J-co 


Let  Xr  =  V;  then 


0(t,y,z)  =  -^dT]  I  dTCp(T,Ti,+0)|2  -  1  dve  ^''^Hi(2)(v) 


.  r2  • 
lJ-00  lJ-00 


(¥) 


With  reference  l8,  page  405,  the  bracketed  term  in  equation  (3^)  can  be 
evaluated  as 


HnV^>'(v)  = 


i— i_L^  1 

r 

t  -  T  ^  , 


MCA  TN  2317 


15 


Thus,  If  the  vortex  of  the  foremost  Mach  cone  Is  taken  as  the  origin, 
there  can  be  written 


0(t,y,z) 


2m'tJ 

Jt 


dr 


(t  -  T)Cp(T,T),-fO) 

\|(  t  -  t)  ^  - 


(38) 


which  is  the  equation  of  the  potential  in  terms  of  the  pressure  coeffi¬ 
cient  on  the  wing.  The  value  of  Cp  is  zero  outside  the  wing.  Since 

the  radical  must  be  positive^  0  is  zero  outside  the  foremost  Mach  cone 
as  predicated  by  other  theories. 


AERODYNMIC  BEHAVIOR  OF  SYMMETRICAL  FLAT  BODIES 
IN  SUPERSONIC  FLOW 


As  the  disturbance  potential  for  the  symmetrical  flat  body  in 
supersonic  flow  has  been  given  in  equation  (l5b)^  with  the  Fourier 
integral  method,  the  present  section  will  show  the  expression  for  the 
pressure  coefficient  and  the  wave  drag  of  various  types  of  wing  plan 
forms  and  airfoil  sections. 


Derivation  of  the  Expression  for  the  Wave  Drag  of  a  Wing 


To  the  order  of  approximation  of  the  linearized  theory,  the  incre- 
.  ment  of  pressure  Ap  anywhere  within  the  disturbance  area  is 


Ap  =  p  -  Po  =  “PUu  =  -PqU 


m^  ht 


(39) 


where  the  relations  x  =  (m2  -  1)1/2^  =  m’t  and  u  =  are  used, 

p  and  p  are  the  local  pressure  and  density,  respectively,  at  the 
point  (x,y,d[:z),  and  p^  and  p^  are  the  free-stream  pressure  and 
dehsity,  respectively . 

To  find  S0o/St,  equation  (l5b)  can  be  differentiated  as 

where  S0o/St  is  an  even  function  with  respect  to  z. 


u 

m* 


bt 
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Substituting  equation  (Ij-O)  into  equation  (39)^ 

Ooo 


Ap(x,y,z)  = 


2m' 


dq 


(41) 


U  -00 


If  the  pressure  increment  on  the  wing  surface  which  is  located 
equivalently  at  z  =  +0  is  desired,  r — >jy  -  q  |  and 

*00 

X  dA,f(X,q)e^4tgQ(2)^^jy  _  ^  (1^2) 

00 


^(t,y,±0)  = 


2m' 


dq 


If  it  is  the  wave  drag  distribution  along  the  span  hB/hj  which  is 
of  interest,  all  the  horizontal  components  of  the  pressure  increment 
along  "tiie  x—  or  "t- direction  on  bolh  "the  lop  and  IdoHohi  surfaces  over  a 
unit  span  width  must  he  summed  up.  Mathematically^ 


ApCq  dx 


(43) 


where  Uq  Is  the  surface  slope  of  the  airfoil  along  the  line  of  con¬ 
stant  y  as  given  in  equation  (18) .  With  the  expression  for  ao(t,y,+0) 
and  ZS)(t,y,+0)  in  equation  (42),  the  drag  distribution  per  unit  span 
can  be  determined  as 


6D  . 


dT\ 

-00  iJ  -00 


X  dXJIo^2^Xly-qi)f(X,q)  f”  dt  f ”  dX' ,y) 


-00  U-00 


(44) 


where  X'  and  y  correspond  to  Oq.  To  evaluate  the  last  double 
integral,  write  X’  =  -X"  and  the  double  integral  can  be  written  as 
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dt  /  ^^f(X'^y)  =  I  dt 


00  /  *00 


(-dX")e^^^‘^"^V-V’,y) 


W  -00 


-00 


-oc  U 


►oo  Ooo 

dt  I  dX"e^(^"^")‘'^f(X'%y) 

-00  U-oo 


where 


f(-^-"^y)  =  fo("^"^y)  +  ifi(-^">y) 


=  fo(^"^y)  -  ifi(>^".y) 


=  f(x",y) 


is  introduced  because  Tq  is  even  and  is  odd.  From  the  complex 

Fourier  integral  theorem, 


2nf(X,y) 


(i^6) 


if  the  Lebesgue  integral 


dx"|f(>t",y) 


<  00  exists  or  if  f(X",y) 


is  of  bounded  variation^  in  the  neighborhood  X  =  X," .  (See  theorem  23^ 
reference  19;  p.  42.)  Substituti^  equation  (46)  into  equation  (44) 
gives • 
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(Dy  ■ 


2flPoU‘ 


00  U”00 


X  dXHo(^)(x.|y  -  n| 


2«PoU‘ 


-00  U-00 


dT]  I  X  dX 


[Ho(2)(>-|y  -  nl)f(^:>n)f(^^y)  + 


Ho^2)(x|y  -  nl)f(>- 


,Ti)f(^.y)] 


=  2JtPou2  I  dT)  I  X  dx/jo()^|y  -  + 

J-oo  Jo  ^ 

fi(x,Ti)fi(5^.y3+ Yo(xly  -  Til)[fi(x,ti)fo(x,y)  -  fo(>^.n)fi(^.y3| 

(1^7) 

where  the  relations 

HQ(2)(_x|y  -  tI)  =  Ho(2)(x|y  -  Til)  =  Jo(My  -  nl)  -  iYo(My  -  tiI) 

f(  -X,r|)f(  -^^y)  =  f(X^Ti)f(X,y) 

are  used  so  that  the  integration  limit  changes  from  0  to  00  instead 
of  -00  to  00.  The  imaginary  part  automatically  cancels  out. 

Now,  if  D,  the  wave  drag  of  the  whole  wing  system,  is  desired, 
equation  (^7)  can  he  integrated  with  respect  to  y. 

p  _  “  2KPqU^  dy  j'  ^  dXHo^^^Xly  -  ■nl)f(X,Ti)f(X,y) 

4 

J-00  J-00  Jo  ^ 

fl(X..Ti)fi(X,y)]  +  Yo(X|y  --  nl )  [fi(X,ii)fo(X,y)  -  fo(X,Ti)fi(A,,  yg} 
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Now^  as  the  wave  drag  D  is  independent  of  y  and  it  should  he 
independent  of  the  order  of  integration  or  of  the  interchange  of  y 
with  Tj.  In  carrying  out  such  an  interchange, 


D’  =  ^irpoU^ 


d\ 


-  y|) 


+ 


+  Yo(x1ti 


y 


(X,y)fQ(X.,T)) 


This  shows  that  the  necessary  condition  D  =  D*  is 


100 

—  00  u 


n«> 


X  dX 


yo(^|y  -  t-iD  [5i(^^y)fo(^^’‘i)  - 


=  0 


This  is  one  of  the  important  relations  that  has  been  developed  by 
Von  Karman. 

There  is  another  way  of  evaluating  the  wave  drag  as  pointed  out  by 
Von  K^rm^n.  The  total  momentum  transfer  through  a  cylindrical  surface 
of  very  large  radius  R  which  contains  the  wing  in  the  flight  direction 
must  be  equal  to  the  wave  drag,  by  means  of  Newton* s  laws.  Thus, 


i2jt 


D  =  -Po  dt 


R  de  MM 
Sr  at 


(50) 


which  is  the  same  as  equation  (48).  Applying  equation  (49),  equation  (48) 
can  be  written  as 


00  nco 


D  =  4nPoU2  dy  I  dr]  X  dXJo(X|y  -  r,  | )  [YqC  X,Ti)fo(  X,y)  + 

ci  -00  U  —00  U  0 

fl(X,ri)fi(X,y^  (51) 
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In  the  case  of  a  sweptback  wing  which  is  symmetrical  with  respect 
to  the  x^z-plane,  fo(k,y)  and  are  even  functions  with 

respect  to  y.  To  take  care  of  this  point,  write 


00  noo 


D  =  4npoU2  dy  dt]  X  dXJo(xly  -  t]  1)  [fo(  X,  1  t)  |)fo(  X,  |y  |) 


-00  U  0 


fl(X,hl)fi(X,|y|)] 

poo  poo  p  00  ^ 

=  knPoU^  J  dy  J  dtj  J  X  dX<|^Jo(X|y  -  t|  | )  + 

Jo[x(y  +  T))]|  [fo(X,Ti)fo(X,  |y  I)  +  fi(  X,r|)fi(  X,  |  y  jfj 

Jnoo  poo  p  00 

dy  I  dri  I  X  dxijo(X|y  -  tiD  + 

0  Jo  J  0  ^ 

Jo[x(y  +  [fo(X,Ti)fo(X,y)  +  fi(X,r,)fi(X,y)]  (52 

where  the  limits  of  rj  and  y  change  as  shown.  Equation  (52)  will  be 
used  in  the  next  section  very  often.  If  the  span  is  finite  (say  that 
b  is  the  semispan)  and  the  airfoil  sections  are  similar  everywhere, 
equation  (52)  can  be  written 


Ob  Ob  n  00  r 

D  =  8jtPoU2  I  dy  dr,  X  dX  j  Jo(  X|  y  -  ti  1 )  + 
I  I  I  A  ^ 


Jol^(y  [/o( x,y)  +  fi( X,Ti)fi(  X,y^ 
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Wave  Drag  of  a  Tapered  Sweptback  Wing 

If  the  double-wedge  profile  is  taken  as  the  airfoil  section  of  the 
wing  as  shown  in  appendix  k,  and  equation  (a8)  is  introduced  into 
equation  (53)^  there  results 

pb  pb  p« 

D  =  SjiPqU^  dy  I  dT)  I  X  dX 

Jo  Jo  Jo 


n  (a 

-  cos  Xay*)(l  -  cos  Xa-q^)cos  Xp(|y|  -  \r\\) 

j  X 


For  details  of  the  notation  refer  to  figure  1  and  appendix  B. 
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8poU‘ 


The  wave  drag  equation  can  he  written  in  another  form 


Ph 


jtD 


Pb' 


dy' 


■(ao')^(ao')^  Jo  Jo  Jo 


Ooo 


dt)^ 


dX 

X 


^<doi^n'  -  y' 


)  ^0*J 


JoO’^’  y‘)^o 


1  -  cos  Xao'[2  +  3(q'  -  y’ ) 


v'(ti'  +  y')J  ^  ^  |l  -  cos  Xao'[2  -  p(q'  -  y'  )  - 

v’(ti'  +  y')J  ^  -  ^|l  -  cos  Xao’CP  +  v’)(il'  -  y' )J  - 

-  cos  [Xa'(p  -  v')(ti'  -  y*)]^  + 

i|i  -  cos  xao’li  +  P(n’  -  y')  -  v'y'3|-  + 

l-|l  -  cos  Xao’[l  -  P(ti’  -  y’)  -  v'y’J^  + 

-  cos  AaQ*[i  +  P(t]’  -  y' )  -  v'Tj^J’ 

i  |l  -  cos  Xao'[l  -  P(q'  -  y' )  -  " 

[1  -  cos  ^Xao'(T1'  -  y'Oj 


(56) 
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With  the  known  infinite  integrals  of  Bessel  functions  as  shown  in 
appendix  C,  the  integration  in  eqtiation  (56)  can  he  carried  out. 


dy-  av  - i  cosh-1  -y">  .  1  cosh-i  -v-(h’  ^y) 

^  Ti*  —  ^  Tl*  + 


-1  g  -  -y*)  -  y*)  _  1  „osh-i  g  -  -  y')  -  v*(n*  y') 


8poUg(ao')2(ao')2  Jo 


i  cosh“^  ^  i  cosh"^  (^  +  V*)  -  V  cosh“^  (p  -  v' )  -  f  cosh 


1  _i  (P  -  v’)(il»  -  yO 

f  cosh  ^  - - - 

4  «»  +  v» 


.1  1  +  P(tj* 

-  y') 

-  v’y' 

1 

cosh-l  ^ 

+  P(ti’  -  y' ) 

-  v*y' 

,  1  ,,.h-i  1  -y)  -'"y 

L 

1’ 

-  y’ 

2 

n’  +  y’ 

2  7)*  +  y' 

.1  1  -  eCh’ 

-  y’) 

"  v'y’ 

+  1 

cosh-1  " 

+  P(ti*  -  y' ) 

-  v*n* 

.  1  cosh-1  "  " 

‘  V 

-  y’ 

2 

n"  -  y’ 

2  V  +  y’ 

1  -  3(ti’ 

-  y’) 

-  V'Tl' 

+  1 

-  P(ti'  -  y’ ) 

-  v'T 

coch-'i  coch--'  B 

1  ' 

V 

-  y' 

^  2 

CUSQ  “ 

r,»  +  y* 

Ti»  +  y» 

In  all  cosh*^  (A)  terms,  A  must  be  positive.  In  order  to  carry  out 
the  integration,  a  transformation  of  the  coordinate  system  may  be  intro¬ 
duced  as  follows:  Let  u’  =  t]’  +  y^,  v’  =  -  y’.  Then, 

S( y*  * Ti’ )  1 

dy^  dT]*  =  — - !—  du'  dv’  =  --  du^  dv^ 

‘  5(uSv»)  2 


y^  =  0,  T)*  =0  — >  u*  =  0,  V*  =  0 

y*  =  0,  Tj*  =  b  — >  u’  =  b,  V*  =  b 

y^  =  b,  Tj*  =0  — >  u*  =  b,  V*  =  -b 

y^  =  b,  T]*  =  b  — >  u*  =  2b^  ^  v  =  0 


With  this  transformation,  the  square  domain  in  the  y *,11* -plane  will 
transform  to  a  diamond-shaped  domain  (of  twice  the  area)  in  the 
u*,v* -plane. 
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_ _JiS_ — =  -(b 

8poU2(a^’)2(ao*)2 


Ob'  pu*  / 

)2ji  cosh-1  (p  +  V')  +  i  co6h-l  (p  -  v’)  +  cosh’l  pj.+  I  du*  I  [  I  P)  - 


/2  +  Pv' 

1  pneb“d 

f2  -  v'u'  _  J 

-  i  cosh-! 

'2  -  .  V') 

U’  ] 

"  ij;  cosn 

k  v'  ^ 

^  U>  / 

X  .,,,.X  (L^],  [.  1  .„3h-x  ^  f)  -  i  CO.K-X  (1^  -  V')  - 


-  ~  cosh-! 

f2-Pv'  -v*) 

-  i  cosh-! 

(p  +  v')v 

4 

\2b'  -  u’  / 

2b’  -  u* 

n  .u  !  Ln=>,-1  ( 

2  -  v'u'  «  .  V'\  r 

2  +  (2p  +  v' )v’  v'l 

-1  +  —  <cosh  1- 

J  2]^  V 

2v’  2/  [ 

2u'  2  J 

/2  -  v'u’  o ...  v*^  ^  „.,.v,-!  r?. 

-  (2p  -  v' 

)v’  v'"] 

+  cosli^^  j 

f2  -  v'u’  p  _V’U 

V  2v’  ^2/  L 

2u’ 

■  2J 

^  2v’  ^2} 

(2  +  (2p  -  V’  )v’  V'1  . 

(2  -  V’u' 

-  P  - 

h-1 

[2  -  (2p  +  v’)v'  v'l 

L  2u’  2J 

V  2v’ 

^  2) 

L  2u’  2J  ■ 

.1  12  -  (2p  -  v')v’  v’l 

/2  -.2v’l3'  +  v'u’  ,  „ 

.ill 

1  [2  +  (2p'  -  v' )v'  v'l 

L  2(2b’  -  u')  2J 

\  2v’ 

2  / 

L  2(2b'  -  u’)  2J 

T  /2  -  2v’'b’  +  v'u' 

-H)  .  cosh-1  P  - 

.vlll 

.  /.cosh-l  Sil  -  cosh-1  f  , 

\  2v» 

2)  ,  L  2(2b’  -  u’  ) 

.Jj 

^  u'  2b'  -  u' 

The  details  of  the  integration  are  too  complicated  to  be  given  here 
The  result  can  be  summarized  in  equation  (59)  which  follows,  where 

Pq  =  p  +  V',  Pi  =  3  -  V,  a  =  1  -  v'b',  and  sin“l  (a)  =  (sign  A.)| 

if  I  A1  >  1.  The  upper  line  with  1  -  3^  should  be  used,  if  3  <  1; 
the  lower  line  with  p^  -  1  should  be  used,  if  3  >  Ij  and  the  signs  i 

of  all  cosh"!  terms  should  be  reversed  if  3  <  -1.  The  above  rules 
should  be  also  applied  for  Pq  and  Pi.  Therefore  equation  (59)  covers 

all  cases  except  a  few  limiting  cases:  3  =  ±1,  Pq  =  ±1;  Pi  =  -1^ 
v'  =  0,  and  p  =  00  (l.e.,  M  — ^  l).  These  limiting  cases  will  be  shown 

in  the  next  section. 


4. 
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L  2{Po®  -  i)^ 


^  r  ^  -  COSll"^  -^  +  1»  COBh'^  —  I  + 


[  A^-  O*’'] 

1.  . 

_j  opo  +  2b’  1 

[;°  2(i  +  5b’jJ 

1  2(00  ♦  5)(l  - 

®  "  a  +  25ob’  "  “  I 

P  <7  *  2Pob’  J  j 

j  1  (1  -  5o^)b’  j 

1  ^  (a  +  25c,b’)  2  1 

r„.h-i  |‘’Po*2b’| 

1 .  L .  {1  *  n 

1  1  °  2(1  ♦  5b’)  1 

J ' .(.. . .)(..» -  1 

L  |o*25ob'! 

1  1  0  +  2Pob'  IJ 

(•  -  f., 

.1  ®®1  ■ 

2(l-pb’)Jj 

2(5i  ♦  5)(l  -  51^)^''^  ]_ 

a  -  25ib’  "  ' 

a  -  26b’  J 

1  L 

1 

!o6i  -  2b’ 1 

1  1  (i-6i®)b’n 

[  1^1  2(1-Pb')| 

J  2(5i  S)(3i2  -  1)^/2 

|_  |e  -  25ib'  j 

1  *  ”■“>  [“i  -  „  .  apb’  IJ 

°-P  •  g*’’  1  .  cosh-1  L 

(1  -  52)  b’ 

(a  +  26b' )2  r  Ib^  +  2b’  ]  ; 

.  (1  -  62)b’  ■ 

0  -  26b' 1  1 

a  -  25b’ 

2(5i  .  P)(52  -  1)1/2  L  1"  "  1 

a  *  25b’ 

- .  Jsln-1  -  Bl^-1  0  -  L  .  - ■  hin-I  SLtJ^  .  ,1„-1  L  ^  (l  -  I 

- - Lsh-l  ISS^L^I  .  cosh-1  L  -  iLjLi!hl|]  .  - f£_LlBbl)f - r  i  ^  .(l  ..pa)b’  l] 

2(00  ♦  0)(02  -  l)l/sL  l^-SPb'l  I  a  -  25b'  IJ  2(Bj  +  P)(b2  ,  1)1/2  L  |b  +  25b’  |  |  b  +  25b’  |J 

__^  [,. .  IL:^] .  ^ ....  [,. .  li^] .  ^ ....  .  (. .  ^  .„.  [^ .  (. .  .... 

;7(;i^  "l)l7^|f  I  -  ""  \^o  -  ~  1  *  -7  l^o  ^  (l  -  1  -  ^  «Bh-l  |5o  *  (l  -  5o2)b’|  .  coBh-1  |5o|j 

[. ....  [.. .  li^] . . ....  [, .  (i^] .  ^ ....  [.. .  (. .  ^ ....  [, . ,. . ^ ....  ,] 

!.-(.iA  i)'*P  I’*  ■  I  *  I*’  •**  I  •  T  I”  ■  t’  -  I'l  -  (‘  -  *i">’'i  ■  I  “"■•*  nJ 

__i^  •■«-’  [•  •  H#i:] .  .,.-1  [, .  t^] .  jfi_ ....  [,.(..  ....  [,,(..  ,a).,] .  ,..-1  ,| 

‘“h-‘  1'  ■  "-  ^-1  •  I’  ■  ^  I'  •  5  ■  ''>•■'1  •  ^  -  «*>'■'  I  -  ,)  I'i] 


where  a,  b,  A,  B,  and  F  designate  the  terms  to  be 

used  in  the  next  section. 

For  a  given  sweptback  wing^  3^  >  or  <  1  means  the  Mach 

line  ahead,  on,  or  behind  the  leading  edge,  respectively;  >  1,  =  1, 

or  <  1  means  the  Mach  line  ahead,  on,  or  behind  the  trailing  edge.  With 
a  fixed  taper  ratio,  the  trailing  edge  may  be  swept  forward;  p^.  =  -1 
means  the  sweptforward  trailing  edge  coincides  with  the  Mach  line. 
Similarly,  p  refers  to  the  midchord  or  the  line  of  the  maximum  thick¬ 
ness.  At  all  P’s  equal  to  |l|  there  occurs  a  discontinuity  of  the 
slope  of  the  curve  of  the  wave  drag.  In  other  words,  ciorves  for  P  =  111 
are  the  envelopes  of  cusps  of  the  curves  at  constant  p. 
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To  demonstrate  the  use  of  these  results,  the  wave  drag  coefficients 
of  two  families  of  sweptback  wings  have  been  calculated  as  shown  in 
figiires  2  and  3  with  taper  ratios  of  0.2  and  0.5,  respectively.  The 
essential  parameters  are  A  tan  A,  and  p.  The  effect  of  the 

Mach  number  is  contained  in  P  =  tan  A/m'  and  Cjjq  =  hcL^jm'  ,  the  two- 

dimensional  wave  drag  coefficient  of  Ackeret.  To  use  these  curves,  first 
fix  the  sweptback  angle  A  and  aspect  ratio.  Of  course,  a  will  have 
to  be  chosen  as  0.2  or  0.5  in  order  to  use  these  graphs.  Then,  with  a 
fixed  A  tan  A,  read  off  at  various  values  of  p  along  the  fixed 

abscissa.  Replotting  C-^jCj^Q  against  p,  the  example  shown  in  figure  4, 
actually  gives  a  family  of  swept  wings  of  fixed  A  tan  A.  The  curve  of  Cj) 

against  Mach  number  can  be  plotted  from  the  relation  M  = 

This  is  not  given  in  this  report. 

It  is  very  interesting  to  note  that  ( /*^Do) max  when  the 

Mach  line  coincides  with  the  leading  edge,  Pq  =  1,  or  the  line  of  maxi¬ 
mum  thickness,  p  =  1.  In  figure  2  (a  =  0.2),  (‘^D/'^Do)max 

p  =  1  for  A  tan  A  >  7-5  approximately,  and  at  Pq  =  1  for 
A  tan  A  <  7.5.  In  figure  3  (o  =  0-5) ^  (Cd/Cdo)  max  always 

at  Po  =  1  except  for  A  tan  A  <  2/3.  Also,  increases  mono- 

tonically  with  A  tan  A  at  Pq  =  1,  P  =  1,  or  P']^  =  l,  but  decreases 
monotonically  with  A  tan  A  for  p^  =  -1. 

Figure  5  shows  p  corresponding  to  Pq  =  1  and  Pj  =  1  plotted 
against  A  tan  A.  This  is  needed  for  locating  the  cusp  of  the  curve 
of  C^yCjjo  against  p,  as  shown  in  figure  h.  These  curves  are  obtained 
from  the  relations: 


2(1  - 
(1  +  a)A 


^) _ 

tan  A 
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A  few  interesting  facts  can  be  drawn  from  figures  2  and  3^  although 
the  taper  ratios  are  fixed  at  0.2  and  0.5^  respectively.  First,  if 
P  =  tan  A/m’  =  0,  M — ^  for  A  /  0,  or  A  =  0  for  M  >  1.  For  A  /  0, 

Cp/Cj3o  =  1  =  ^Do  indicated.  For  A  =  0,  0  <  0j)^Cj)q  =  1 

and  the  exact  value  of  is  indeterminate.  As  p  increases 

from  0  to  1,  rises  quickly  with  A  tan  A  until  it  meets  the 

cuDTves  Pq  =  1  and  then  varies  with  A  tan  A  very  slowly.  At  p  =  1, 
C])|Cj)o  is  extremely  high.  This  means  that  the  closer  the  Mach  line  to 

the  maximum-thickness  line,  the  higher  the  drag.  As  p  increases 
further,  from  1-^  00^  Cd/Cj)o  decreases  with  increasing  p.  This  means 

that  the  closer  the  Mach  line  to  the  maximum-thickness  line,  the  higher 
the  drag  also. 

If  the  sweptback  angle  is  negative,  or  if  the  flow  direction  is 

reversed,  the  curve. of  against  -A  tan  |  a|  is  actually  the 

reflection  of  the  present  curves  with  the  following  changes:  p  — >  -p, 

Po  =  1  — ^  Pi  =  -1,  Pi  =  1  — ^  Po  =  and  p3_  =  -1 — >  Po  =  1  if  Po 
is  defined  as  corresponding  to  the  leading  edge,  and  P]_  as  corre¬ 
sponding  to  the  trailing  edge.  It  will  be  found  that  depends 

only  on  A  tan  |A|  and  lP|  for  the  same  wing  whether  it  flies  forward 
or  backward.  This  confirms  the  reversed-flow  theorem  in  references  2 
and  5* 


Special  Cases  of  the  Wave  Drag  Equation 

Most  of  the  special  cases  of  the  wave  drag  equation  are  required  in 
calculating  the  curves  shown  in  figures  2  and  3  and  therefore  must  be 
written  down  one  by  one.  These  special  cases  are  based  upon  the  taper 
ratio  that  is  between  0  and  1. 

(a)  As  Po  — >1,  expressions  A  and  D  in  equation  (58)  take  on  the 
following  limiting  values  while  the  rest  of  the  equation  stays  the  same 
except  a  =  b  =  c  =  0. 
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(1  -  v't’  {1  -  v'b'  -  b'\l/2  /  1  -  V't'  \l/2 


2(v'  -  1) 


1  -  v'b’  +  b' 


1  -  v'b'  +  b' 


(1  -  v'b'  +  2b'  /l  -  v'b'  -  2b' W2 

2(v'  -  2)  U  -  v'b'  +  2b'/ 


-v'b'  +  2b 


1  -  v'b 


1  -  v'b'  +  2b' 


=  —  -  v'b'  +  b' 

2v'  I  L 


-  (1  -  v'b'  +  2b' 


2V'  -  1  (1  _  ■b.)l/2  _  2(V’  -  1) 


(1  -  2b'  + 


v'  -  1 


V'  -  2 


(V  -  l)(v'  -  2) 


This  case  occurs  when  the  right  characteristic  coincides  with  the 
leading  edge  on  the  right  side. 

(b)  Similarly,  for  =  1,  expressions  B  and  E  become 


(1  -  v'b'  -  b')^  /l  -  v'b'  +  b 


2(v'  +  1) 


1  -  v'b'  -  b* 


>\l/2  /  1  _  v'b'  \^/2 

'}  \l  -  V'b'  -  b'/ 


(1  ■-  v'b'  -  2b'  /l  -  v'b'  +  2b'V/2 

2(v'  +2)  ll  -  V'b'  -  2b'/ 


1  -  V'b 


1  -  v'b'  -  2b' 


=  -  v'b'  -  b')V2  +  (1  -  v'b'  -  2b' _ 


1  (1  +  V)l/2  +  Mll-Lll  (1  +  2V)^/2  +  - - - — 

V*  +  1  v’  +  2  (v^+  l)(v’  +  2) 


Here  b  =  0.  All  the  other  tenus  remain  the  same  as  shown  in  equa~ 
tion  (59) •  This  case  occurs  when  the  right  characteristic  coincides 
with  the  trailing  edge  on  the  right  side. 
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(c)  When  3  =  1, .the  expressions  C  and  F  are 


(  a  +  2b' 

2(2  -  v’ ) 

(g  -  2b' 

2(2  +  v'  ) 


M62) 


(1  -  2b'  )V2  *  V)  (1  ,  2b'  )l/2  - 

2-v'  '  2  +  v'  (1..V.2) 


Here  b  =  c  =  0,  All  the  rest  remain  the  same  as  shown  in  equation  (59)* 


(d)  When  Pq  =  ^  ^  become 


A  =  B  in  equation  (6l),  where  =  1 


D  =  E  in  equation  (6l) 


(63) 


Here  a  =  0.  The  rest  of  the  terms  in  equation  (59)  remain  the  same. 
This  case  occurs  when  the  left  characteristic  coincides  with  the  leading 
edge  on  the  right  side.  This  conforms  with  the  reversed-flow  theorem, 

Pq  =  -1 — ^  3^  =  if  the  flow  direction  is  reversed. 

(e)  When  Pj^  =  -1,  B  and  E  become 


B  =  A  in  equation  (60),  where 
E  =  D  in  equation  (6o) 


(64) 


Here  a  =  b  =  c  =  0.  The  rest  of  the  terms  remain  the  same  as  given  in 
equation  (59)* 
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(f)  When  p  =  -1,  C  and  F  are  the  same  as  given  above  for 
3=1  with  some  rearrangement  of  terms  and  a  =  c  =  0. 

(g)  When  v’  =  0,  or  there  is  no  taper^  3q  =  Pp  =  3  and  equa¬ 
tion  (59)  becomes 


rtf  =  (b«)2 


- — 7~  cosh-1  I - 


^1  -  cosh-1 

+  J)-  COEh“l  l-^j 

1  2P  1  lb’ 

[2b’  1 

-sin-I  .  sln-l 

^  Lh-1  -  co=h- 


I  .  L„-l  .  sln-1  L  . 

L  2(l  +.pV)J|  2p(l  -  p2)l/2  I  1  +  2P1)’  L  1  +  2pt'  J 

V  n  ^  |£_Liv.| .  .c.h-1  L .  CLiifki] 

)  U  2p(p2  -  l)V2  [  11  +  2pb*  1  I  1  +  2pf  IJ 


(1  -  P^)^' 

2(1  +  pb' 


(1  -  Pt’)' 


2p(l  -  p2)^/^ 
(1  -  Pt*)2 


.  1  P  -  "b'  .  1  L  (1  -  P^)^’ 

J-sin~^  -  +  sin“^  P - - r 

1  -  pb*  L  2(1  -PV) 


(1  -  2pb')‘ 


Icosh-l  -  cosh-1  p 


I  »  cQsn^-^  — - 

2p(p2  .  1)1/2  [  |l  -  @1 


2(1  -  Pb') 


_ ^  Ln-1  .  slc-1  I  .  kz^ 

.  y^.Lsh-l  li^l  .  cosh-1  L  - 

2p(p2  -  l)V2  [  il-2pb*|  I  l  +  2pb>  IJ 


2p(l  -  p2)^/2 


J  Bin'l 


2p(p2  -  1)^/2  L 


(1  - 


jp  +  -  sin-1  jp  +  {1  -  p2)bjj  +  sin-1  |p  -  If  "  "  P^)l>'lj 

+ 

-cosh"l  Ip  +  ii— -~f-j  +  cosh'l  |p  +  (1  -  p^)b’  I  -  cosh-l  jp  -  ^  cosh"!  |p  -  (l  -  p2)b’|J 

.  ,1.-1  [3  .  *  2  sln-l  [3  .  (1  -  32)"/^b]  -  2  sih-1  [3  -  (1  -  32)"^^b]j 


sin-1 


- . y-  -cosh’^  |p  -  ~  I  +  cosh'l  L  +  (1  -  P^)p  I  _  2  cosh"^  |p  +  (1  -  p2)^^%'|  -e  2  cosh'l  L  _  (i  _  p2)^/^\ 

(p2  _  1)1/2  |_  I  2  I  I  2  I  I  I  I 


(l  -  ‘ 

(p2  11^ 


r  ^  7^'  ■ 

r  ^  -11/.^ 

1 

1/2 

_(2  +  pb')^  -  (b’)^ 

_(2  -  Pb’)2  -  (b’)2j 

(1  -H  pb*  )2  -  (b*)2 

(1  -  pb*)2  -  (b’)2^ 

4  1  r  4  _  r  1  _  r  i  f^] 

+  Pb*  f  -  (b*  )^J  _(2  -  Pb*  )2  -  (b*  )2j  [(1  +  Pb*  f  -  (b*  )^J  L(1  -  P^’  “  (^’  fj  J 


(65) 


The  above  equation  is  in  accord  with  the  results  given  in  reference  2. 
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(h)  If  M  — >  1^  (3  — ^  ±00  as  A  is  positive  or  negative.  By 
definition^  p.  =  t/ao^  v*b*  =  vp^  o)  =  p  tan  =  p  tan  Aq,  and 

0)^  =  p  tan  A^^. 


There  are  a  number  of  other  limiting  cases  such  as  Pq  +  0  =  0 
(Aq  =  -A)  and  P3_  +  P  =  0  (A]_  =  -A)  which  are  also  required  in  cal¬ 
culation.  If  a  =  0  and  P]_  =  0^  equation  (59)  will  give  the  wave 
drag  of  the  triangular  wing  which  checks  with  Puckett*  s  results.  Owing 
to  limited  space ^  they  are  all  omitted. 


(66) 
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AERODYMMIC  BEHAVIOR  OF  A  WING  WITH  A 
GIVEN  LIFT  DISTRIBUTION 


From  the  section  Supersonic  Flow  about  a  Lifting  Surface, 


g(^^y)  =  ^ 

2jt 


dt  e-^^^Cp(t,y,+0) 


(30) 


Thus,  when  the  pressure-coefficient  distribution  Cp(t,y,+0)  is  given 
on  the  wing  plan  form,  g(^^y)  can  be  obtained  with  this  equation. 

Consider  a  tapered  sweptback  wing  with  a  constant  lift  distribution 


Cp( t ,y)  —  CpQ 

\ 

-b  ^  y  =  b 

Pobl  -  ^o' =  ^  =  Pilyl^o’ 

where  Pq  and  are  defined  in  figure  1.  Thus,  at  any  y  within 

the  span. 


g(  >-,y) 


■^Pi|y|ao' 

_1. 

2jt 

Jpol  yl-ao’ 


-iXt 

Cpoe  dt 


^^po  r -iHPl|y|ao')  ^-iH3o|y|-ao’) 

2jtX  Lf  "  ® 


(67) 


NACA  TN  2317 


33 


As  A.  — >  0,  g(0,y)  can  be  evaluated  as 


g(0,y)  =  — [ao 

JT 


Po  -  Pi 


|y| 


(68) 


Substituting  the  above  equation  in  equation  (24),  there  results 

Pb 


w(t,y,+0)  = 


-m'U 


dq 


-b 


(y  -  Tl)^ 


,  Po  ■  Pi,  , 
n  \^o  ■  2  hij 


|y  -  Til 


Hoo 


dX 


^PO  r-iX(pi|y|+ao’) 

2nXl  l! 


-m'UC 


po 


2jt 


Pb 


-b 


dTi 


(y  -  Tl) 


^ao’  - 


Po  -  Pi, 


i^(Pl|y|+ao’) 


U-  00 


(69) 
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With  the  aid  of  appendix  the  infinite  integrals  can  be  evaluated 

and 


ly  -  ^1 
a 


+  0,  ti  -  3ihl  ^  ly  -  nl 

|y  -  tl 


ti  -  Pihi  ^  -  -  (y  - 

|y  -  Til  jy  -  ’ll 


Pilnl  &  ly  -  'll 


where  to  =  t  +  aoS  **"1  •  Only  one  of  the  two  expressions 

holds  for  the  range  specified  in  the  bracket. 


MCA  TN  2317 


If  the  integration  of  is  written  from  0  to  h  only, 

[^^•o'  “  (Po  "  Pl)3  23^'  -  -  P2.)'n 

(y  -  Tl)^  (y  +  tj)2 

gap'  -  (Po  -  Pi)ti  ^  2ao'  -  (Po  -  Pi)^ 

(y  -  3)^  (y  +  r,)2 


2\/(^o  -  Po^)^  -  (y  -  ti)2  2^ 

f^o  -  Po'^)^  -  (y  + 

OJ 

l>i 

(y  +  n)2 

2  ]l(ti  -  Pin)^  -  (y  -  ti)2  2  \j 

(ti  -  PiTi)2  _  (y  +  ti)2 

ro 

(y  +  n)^ 

w(t,y,+0)  = 


II 


to  -  Pon)^  -  (y  -  n)^ 

(y  -  n)^ 

\|(^i  -  -  (y  -  n)^  ^  \/ ("^o  -  Pq^i)^  -  (y  +  3)^ 

(y  -  3)^  (y  +  t])2 


Ph  r 


•Cpo^’U 

2jt 


JO 


ii 


\[(^i  -  -  (y  + 

(y  + 


where  to  -  =  2ao*  is  used  so  that  the  first  terms  cancel  out. 
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Now,  carrying  out  the  integration  with  respect  to  t). 


w(t,y,+0) 


Cpo®'  U 

2it 


[-(A  -  B)  +  (C  - 


(72) 


where 


rib 


A  = 


dTl 


Nj(to  -  Pon)^  -  (y  - 
(y  - 


(see  table  l) 


(73) 


B  = 


ph 


Uo 


dri 


\[(tp  -  -  (y  - 

(y  -b)2 


(see  footnote,  table  I) 


(7i^) 


nb 


c  = 


-dti 


Jo 


^(to  -  Po^)^  -  (y  + 
(y  + 


(see  table  II )  (75) 


D  = 


_ 

\  (^1  ■  -  (y  + 

JL — -  (see  footnote,  table  II) 


-dT] 


Uo 


(y  + 


(76) 


The  present  treatment  is  equivalent  to  considering  two  hypothetical 
wings.  The  first  one  is  of  constant  positive  lift  distribution  (Cpo  is 
an  assigned  constant)  starting  at  the  leading  edge  and  extending  down¬ 
stream  to  infinity  with  semispan  b.  The  second  wing  is  of  constant 
negative  lift  distribution  (Cpo  is  an  assigned  negative  constant) 
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starting  at  the  trailing  edge  and  also  extending  downstream  to  infinity 
with  semispan  b.  The  superposition  of  both  hypothetical  wings  gives 
the  aerodynamical  behavior  of  the  actual  wing. 

Integral  k,  equation  (73)^  gives  the  contribution  of  downwash  due 
to  the  right-side  area  of  the  first  hypothetical  wing  in  the  forward 
Mach  cone  of  the  point  P(t,y).  Similarly,  integral  C,  equation  (75)^ 
gives  the  contribution  for  the  area  on  the  left  side. 

On  the  other  hand,  integral  B,  equation  (7^),  gives  the  contri¬ 
bution  of  the  right  side  of  the  second  hypothetical  wing.  Similarly, 
Integral  D,  equation  (76),  gives  the  left-side  contribution. 


As  far  as  the  limits  of  the  integrals  are  concerned,  each  holds 
for  within  the  limit  0  to  b,  whenever  the  integrand  is  a  positive 
real  quantity.  However,  for  many  locations  of  P(t,y),  the  integrand 
is  positive  only  in  a  much  narrower  range  of  q  than  the  interval  (0,b). 
In  order  to  determine  the  valid  range  of  the  integrals  and  their  respec¬ 
tive  values  for  all  possible  locations  of  P,  tables  I  and  II  are  given. 
The  corresponding  values  of  the  integrals  for  the  different  cases  are 
given  in  the  right-hand  coliunns  of  the  tables. 

Take  table  I  as  an  example.  If  the  leading  edge  is  subsonic 
( Po  >  l)  ^  'the  lower  limit  of  Tj  is  always  zero  because  the  forward 

Mach  cone  at  P  always  Intersects  the  wing  center  line  within  the  wing 
area.  When  Po  >  1  the  upper  limit  of  t]  is  given  in  the  upper  half 
of  table  I.  Condition  1  concerns  the  location  of  the  point  P,  which  may 
be  ahead  of,  on,  or  behind  the  leading  edge,  while  condition  2  concerns 
the  right  intersection  of  the  forward  Mach  cone  of  P  with  the  leading 
edge,  or  side  edge.  Under  both  conditions  1  and  2,  there  are  six 
possible  cases.  The  upper  limit  of  q  is  given  in  each  case,  and  the 
value  of  the  integral  A  is  also  given  accordingly. 


If  the  leading  edge  is  supersonic  (Pq  <  l)^  the  integral  A  exists 
only  when  the  point  P  is  behind  the  leading  edge.  In  particular, 
both  the  upper  and  the  lower  limit  may  vary  within  the  interval  (0,b). 
The  condition  for  the  lower  limit  depends  on  the  Intersection  of  the 
left  side  of  the  forward  Mach  cone  of  P.  It  is  zero  when  and  only 
when  the  intersection  is  on  the  center  line  of  the  wing  in  the  wing 
area.  If  the  intersection  is  on  the  right-side  leading  edge,  the  lower 
y  -  tQ 

limit  is  - - —  as  indicated.  The  condition  for  the  upper  limit 

^  “  Po 

depends  on  the  right-hand  intersection  of  the  forward  Mach  cone  of  P. 
If  the  right  side  of  the  forward  Mach  cone  cuts  the  leading  edge  within 


the  span,  the  upper  limit  is 


tp  +  y 

1  +  Po 


;  otherwise,  the  upper  limit  of  r\ 
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The  tables  for  the  other  three  integrals  can  be  explained  similarly. 
Figure  6  shows  four  typical  cases  concerning  the  integral  limits.  Fig¬ 
ure  6(a)  shows  the  forward  Mach  cone  of  a  point  P(t,y)  that  lies  only 
in  the  right-side  area  of  the  wing.  Therefore  only  integral  A  exists 
and  all  the  other  three  integrals  are  zero.  Since  the  leading  edge  is 


/  \  to  +  y 

supersonic  Oo  <  the  upper  integral  limit  is  - .  It  belongs 

^  1  + 
to  integral  A-9  as  indicated,  according  to  table  I.  Figure  6(b)  shows 
the  forward  Mach  cone  covering  both  the  left-  and  right-side  areas  of 
the  wing.  Thus,  only  integrals  A  and  C  exist.  The  integrals  are  indi¬ 
cated.  Figures  6(c)  and  6(d)  can  be  similarly  explained. 


Although  the  lift  distribution  of  the  wing  with  supersonic  trailing 
edge  is  well-»-known,  the  downwash  may  be  interesting  to  explore.  Fig¬ 
ure  7  shows  three  infinite  half -wings  with  trailing  edges  at  15*^^  30°, 
and  45°  from  the  leading  edge  which  is  normal  to  the  direction  of  flight. 
This  is  shown  in  the  t,y-plane,  or  the  x,y-plane  with  a  Mach  number  equal 
to  A  negative  infinite  downwash  always  occurs  at  the  tip.  The 

curves  are  plotted  as  27tw(y/t)ym^CpQU  against  y/t.  As  is  known, 

CpQ  =  2a2_/m’ ,  where  is  the  angle  of  attack  of  the  wing;  therefore, 

jrw/a2_U  can  be  plotted  against  y/t.  Between  the  tip  cone  and  the  wing 

the  downwash  is  constant  but  increases  with  increasing  trailing-edge 
angle.  Owing  to  the  conical  flow,  the  downwash  is  identically  the  same 
along  the  radial  lines  expressed  in  the  conical  coordinate  y/t.  The 
above  curves  were  calculated  with  equation  (72)  and  the  tables  for 
integrals  A  and  B  by  setting  Pq  =  0;  3i  =  tan  15°,  tan  30°,  and  tan  ^5^; 
and  a^*  =  0.  (Here  t  =  to  =  t^- ) 

Figure  8  shows  a  wing  tip  of  unit  chord  with  a  raked  angle  of  30° 
for  three  values  of  t.  At  t  =  0.5,  it  behaves  exactly  the  same  as 
the  30°  case  of  figure  7(b),  but  at  t  =  1.5,  it  is  quite  different.' 

Two  infinite  downwashes  occur,  one  at  the  center  of  the  leading-edge 
tip  cone,  and  the  other  at  the  center  of  the  trailing-edge  tip  cone. 

The  former  remains  negative,  the  latter  positive.  Outside  the  two  tip 
Mach  cones,  the  downwash  is  identically  zero  as  predicted  by  the  two- 
dimensional  theory.  The  above  calculations  can  be  made  from  equation  (72) 
and  the  tables  for  integrals  A,  B,  and  C  by  letting  Pq  Pi  =  30° 

and  ao*  =  0.5*  (Here  t  =  to  -  0.5  =  tj^  +  0.5.) 


As  another  interesting  example,  the  downwash  distribution  near  the 
nose  of  an  infinite  constant-chord  sweptback  wing  is  calculated  (fig.  9) • 
The  leading  edge  is  swept  back  at  30°  but  is  still  supersonic.  The 
downwash  at  different  locations  on  the  span  is  plotted  along  t.  In 
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this  case^  the  curves  are  plotted  with 


at  different  values  of  y  where  CpQ  = 


At  y  =  0,  the  downwash  angle  reaches  the  maximum  value  (O.8165) 
at  a  single  nose  point  and  remains  zero  along  the  chord  up  to  the 
trailing  edge.  There  the  downwash  angle  is  negative  and  equal  to  the 
maximum  value  in  magnitude  and  recovers  to  zero  downwash  downstream. 

For  instance,  at  y  =  0.4,  the  downwash  angle  is  constant  in  the  super¬ 
sonic  region,  and  drops  down  under  the  influence  of  the  nose  Mach  cone. 
At  the  trailing  edge,  the  downwash  angle  drops  to  a  negative  value 
abruptly  and  continues  to  decrease  until  the  trailing  Mach  cone  is 
reached.  The  drop  of  downwash  at  the  trailing  edge  is  exactly  equal  to 
the  downwash  angle  at  the  supersonic  leading  edge.  Further  along  t, 
the  downwash  angle  rises  again  and  becomes  asymptotic  to  0  as  t  — >  ». 
The  same  can  be  applied  at  any  value  of  y  until  the  nose  Mach  cone  is 
off  the  trailing  edge. 

Of  course,  it  is  very  difficult  to  build  a  wing  with  the  angle  of 
attack  as  required  by  the  constant  pressure  distribution.  But  by  the 
present  method  it  is  so  simple  to  calculate  the  downwash  distribution 
anywhere  in  the  plane  of  the  wing  that  it  will  surely  afford  many 
applications  if  the  trailing  edge  is  supersonic.  With  a  given  pressure 
distribution,  it  is  easy  to  calculate  the  sidewash;  however,  the  details 
are  omitted  here. 

As  another  example,  the  downwash  of  a  finite  triangular  wing  with 
a  60°  sweptback  angle  is  calculated.  Figure  10  shows  the  downwash  dis¬ 
tributions  at  different  values  of  y.  Figure  11  shows  the  downwash 
distribution  of  a  tapered  sweptback  wing.  The  taper  ratio  is  1/6. 

Figure  12  shows  the  downwash  of  the  same  sweptback  wing  far  downstream. 

It  is  very  interesting  to  examine  the  effect  of  different  Mach  cones 
on  the  downwash.  One  important  feature  of  the  subsonic  leading  edge  is 
the  infinite  downwash  angle.  As  before,  the  above  curves  are  calculated 
with  equation  (72)  and  the  tables  for  integrals  A,  B,  C,  and  D. 

The  present  approach  to  calculating  downwash  of  a  wing  with  constant 
lift  distribution  cannot  be  applied  to  wings  with  subsonic  trailing 
edges,  because  the  Joukowskl-Kutta  condition  must  be  satisfied  at  such 
trailing  edges.  In  order  to  investigate  the  case  with  a  subsonic 
trailing  edge,  it  is  necessary  to  assume  the  distribution  of  the  pressure 
coefficient  drops  to  zero  at  the  trailing  edge. 


~  v(t,y) 


2a- 


U 


against  t 


'\fl  -  Po' 


is  introduced. 


m 


ko 
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The  preceding  discussion  just  gives  some  cases  with  simple  pres¬ 
sure  distribution  along  the  chord  to  demonstrate  the  application  of 
the  method.  Of  course  the  method  is  not  limited  to  these  cases  alone. 


The  Johns  Hopkins  University 

Baltimore,  Md. ,  November  4,  19^9 
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APPENDIX  A 


FOURIER  TRANSFORM  OF  THE  SLOPE  DISTRIBUTION  ON  AN  AIRFOIL 

As  an  example  of  the  slope  distribution  on  an  airfoil,  the  Fourier 
transform  of  the  slope  on  a  douhle-wedge  airfoil  is  calculated. 


z 


The  accompanying  figure  shows  the  shape  and  slope  of  a  douhle-wedge 
airfoil  on  a  swepthack  wing  at  a  distance  y  from  the  center  of  the 
wing  span.  The  slope  distribution  is  as  follows: 


a-(t,y)  =  0,  t  <to(y)  -  a'(y) 
a(t,y)  =  ao’(y),  to(y)  -  a'(y)  <t  ^to 
a-(t,y)  =  -oto'(y)j  to  <  t  <tQ(y)  +  a'(y) 
a(t,y)  =  0,  t  >to(y)  +  a'(y) 


(Al) 
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For  a  rectangular  wing  with  constant  airfoil  section,  where  a^' , 
Iq,  and  a’  are  independent  of  y  in  the  span,  equations  (a4)  and  (A5) 
reduce  to 


=  0 

a_’ 

=  -r-  (1  -  cos  Xa') 
jtX, 


if  to  is  chosen  equal  to  zero. 

For  the  case  of  a  straight  tapered  sweptback  wing,  as  shown  in 
figure  1, 


a'  (y) 


ao  -■  v|y| 


% 

■to(y)  =  3|y| 


0^0’ (y)  =  independent  of  y 


For  the  details  of  the  definitions  of  the  different  parameters^  refer 
to  the  figure.  Then, 


fo(^^y)  =  —  sin  (X,3|y(){l  -  COS  Xay*) 
cos  (Xp|y|)(l  -  cos  Xay’) 


> 


(A6) 
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APPENDIX  B 

SIMBOLS 

a 

velocity  of  sound 

^0 

half  root  chord 

®o’ 

=  ao/“i' 

half  tip  chord 

^t' 

=  at/m’ 

half  chord  at  y  ^a^  -  v ly|) 

"y’ 

=  ay/m' 

A 

(  2b  \ 
aspect  ratio  ( - ) 

V^o  +  H} 

"b 

half  span 

^  A  tan  A 

ao*  2p 


Cd 

^Do 

Cp 

^po 


vave  drag  coef fig lent 

two-dimensional  wave  drag  coefficient  m’) 

lift  coefficient 
pressure  coefficient 

two-dimensional  pressure  coefficient  or  assigned  pressure 
coefficient 


46 

F(  X,,Ti) 

f(X,Tl)  = 

G(  ) 

g(M)  = 

) 

) 

) 

Hi(2)(  ) 

Hi(2)(  ) 

Jo(  ) 

Jl(  ) 

M 
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distribution  function  of  the  harmonic  source  (complex); 
see  equation  (5) 

-  I  F(M) 

doublet  distribution  function  (complex);  see  equation  (7) 

±XG{  Xyr\ ) 

8m* 

Hankel  function  of  the  first  kind  of  zero  order 

Hankel  function  of  the  second  kind  of  zero  order 

Hankel  function  of  the  first  kind  of  first  order 

Hankel  function  of  the  second  kind  of  first  order 

conjugate  of  H3_^^^(  )  (j]_(  )  +  ±Yi{  )^ 

Bessel  function  of  zero  order 
Bessel  function  of  first  order 
Mach  number 


m'  =  -  1 

n  unit  vector  normal  to  solid  surface 

p  pressure 

Pq  free-stream  static  pressure 


q 


velocity  vector 
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r 


^(y  -  n)^  + 


t  longitudinal  coordinate  in  transformed  plane  (equivalent 

to  time  in  supersonic  case)  (x/m* ) 

to  =  t  +  Bq’ 
tp  =  t  -  ao* 


u 

free-stream  velocity 

u 

velocity  component  in 

X-  (or  t-)  direction 

V 

velocity  component  in 

y-direction 

w 

velocity  component  in 

z-direction 

X 

longitudinal  axis 

y 

spanwise  axis 

y'  =  y/so' 
Yo(  ) 

Yi(  ) 

Z 

a 


a. 


a- 


Bessel  function  of  the  second  kind  of  zero  order 
Bessel  function  of  the  second  kind  of  first  order 
vertical  axis 
angle  of  attack 

angle  of  attack  of  upper  surface  of  airfoil 
angle  of  attack  of  lower  surface  of  airfoil 
slope  of  upper  surface  on  symmetrical  airfoil 
angle  of  attack  of  cambered  wing 


tan  A 


tan  A’ 


Po  ==  P  V* 
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APPENDIX  C 


INTEGRALS  FOR  USE  IN  TEXT 


The  following  integrals  are  those  which  are  used  often  in  the  text 
and  many  of  which  are  not  available  in  the  ordinary  handbook.  In  the 
first  13  formulas^  the  upper  line  corresponds  to  3  <  1,  and  the  lower 
line  to  3  >  1, 
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^  Yo(®t)(l  -  cos  -bt)  = 


K“-*  i1 


(a  ^  "b) 


(b^a) 


y2(at)  cos  bt 


(a  >  b) 


(b  >  a) 


Yi(at)  sin  bt 


1  J>.2  _  „2  _  2 


(a  >  b) 


(b  >  a) 


pa 

Jo 


t«2s+2 
•  a  jf 


^  Y3^(at)(l  -  cos  bt) 


(a  >  b) 


V(ti  -  h^f  -  (y  -  n)2  ^  ti  -  3in 

dj)  -I - - - =  djl  =  ^  +  p,  cosh”l  -^= - —  + 

(y-Ti)2  v2  V  y-n 


2 —  1  (^1^  -  1)11  +  y  -  ^ltl 

,2-1  cosh-l  i - — 

^1  ~  PlY 


1  -  31^  sin-l 


(Pi^  -  i)n  +  y  -  Piti 


ti  -  Piy 


(y  +  n)' 


(ti  -  -  (y  +  n)^  ^  ^  ti  -  PiTi 

.  .  .  .,  =  —  dn  =  -  -3— —  -  6,  ■'•  '■  ■  . . 


J__  dJi  =  -  J- - p,  cosh"-*-  - 

V » 2  V  y  +  Ti 


■ti  +  ^ly 


-Ml  -  sin 


.1  (Pi^  -  i)i  -  (y  H-  Pj,tJ 


ti  +  Piy 


In  the  last  two  formulas,  use  upper  line  if  Pp  >  1,  lower  line  if  Pp  <  1 
Similar  integrals  for  Po  are  obtained  by  substituting  to  for  tp, 

Po  for  pp. 
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TABLE  I 

DPTEGHATION  LIMITS  OF  i)  AKD  CORBESPONDDIO  VALUES  OF  THE  IMTEGRAL  A^ 
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UnSGRATION  LIMITS  OF  AMD  COBBESPONDING  VALUES  OF  THE  INTEGBAL  C'* 


Po  >  1  CThe  Integral  exists  only  when  to  >  y;  the  lower  limit  is. always  zero.) 


Condition  Upper 
limit 


2  b 


+  3^  eosh-1  ^  cosh-l  - 

y  y  K  O  I  to  +  Poy 

/(to  -  -  (y  .  b)^  _  (^o^-i)r..(y.,,t^ 


^o  +  Poy 


Po  <1  (The  Integral  exists  only  when  t  >  y;  the  lower  limit  is  always  zero.) 


Condition  Upper 
limit 


b-:.  .  3,  oooh-i  .  JTr^  -  p°^)  -  ^  goto) 

y  +  t  y+bN  to  +  ^oy 
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(a)  In  physical  plane  (x,y). 

Figure  1.-  Sweptback  wing  in  physical  and  transformed  planes 


PARAMETER  RELATIONS  BETWEEN  PHYSICAL  PLANE  (x,y) 
AND  TRANSFORMED  PLANE  (t,y) 


(b)  In  transformed  plane  (t,y).  Mach  line  at  45 
Figure  1.-  Concluded. 
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Figure  2.-  Wave  drag  coefficients  for  the  family  of  sweptback  wings  with  taper 

ratio  of  0.2. 


NicO 


Go 


Figure  3.- 


Wave  drag  coefficients  for  the  family 

ratio  of  0.5. 


of  sweptback  wings  with  tap 


Location  of  cusps  =  1  and  P-j^  =  of  the  wave  drag  parameter  Cq^Cj-jq  in  terms 

and  A  tan  A  at  a  =  0.2  and  0.5. 
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Trailing  edge  at  15°  from  leading  edge  which  is  normal  to  direction 

of  flight. 

Figure  7.-  Downwash  distribution  for  three  infinite  half -wings. 
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-2  -1  0  1  2  3  4  5 


y/t 

(b)  Trailing  edge  at  30°  from  leading  edge  which  is  normal  to  direction 

of  flight. 


Figure  7.-  Continued. 


(c)  Trailing  edge  at  45*^  from  leading  edge  which  is  normal  to  direction 

of  flight. 

Figure  7.-  Concluded. 


nw 
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Figure  8.-  Downwash  distribution  over  a  wing  tip  of  unit  chord  with  a  raked 

angle  of  30°  for  three  values  of  t. 
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Figure  11.-  Downwash  distribution  at  different  values  of  y  of  a  tapered  sweptback  wing 
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Figure  12. -  Downwash  distribution  at  different  values  of  y  of  a  tapered  sweptback  wing  far  downstream. 


